Sterile neutrinos, dark matter, and the pulsar velocities in models with a Higgs singlet 
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We identify the range of parameters for which the sterile neutrinos can simultaneously explain the 
cosmological dark matter and the observed velocities of pulsars. To satisfy all cosmological bounds, 
the relic sterile neutrinos must be produced sufficiently cold. This is possible in a class of models 
with a gauge-singlet Higgs boson coupled to the neutrinos. Sterile dark matter can be detected by 
the x-ray telescopes. The presence of the singlet in the Higgs sector can be tested at the Large 
Hadron Collider. 
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Sterile neutrinos with masses of the order of a few keV 
can explain the observed velocities of pulsars [l], 0, 0, 0] , 
can play a role in the star formation and reionization 
of the universe Jm , and can be the cosmological dark 
matter 0, 0, 0]if their relic population is sufficiently 
"cold" . 

While it is possible that the keV mass scale is a new 
fundamental constant of nature, it is of interest to con- 
sider the possibility that the sterile neutrino masses arise 
from the Higgs mechanism, just like the masses of the 
other fermions. Then the keV scale is generated by the 
Higgs vacuum expectation value (VEV), and it depends 
on the coupling. The sterile neutrinos can couple to 
the SU3XSU2XU1 singlet Higgs boson S, whose vacuum 
expectation value can give them the Majorana masses. 
Models of this kind have been proposed [13, EH, an d 
their implications for the electroweak phase transition, 
baryogenesis, and collider searches have been studied in 
detail 12J. The production of relic sterile neutrinos has 
also been studied in one specific limit of parameters, in 
which the singlet S with a sub-GeV mass can play the 
role of the inflaton [ll[ . Here we will consider a more gen- 
eral range of parameters, focusing on the region in which 
the sterile neutrinos can simultaneously explain pulsar 
kicks and dark matter, while satisfying the Lyman-a, x- 
ray, and other bounds. We will see that, if the scalar 
Higgs mass is of the order of the electroweak scale, all 
these constraints can be satisfied simultaneously, and, in 
particular, the momentum distribution of the relic ster- 
ile neutrinos can be cold enough for dark matter. The 
mixing between the Higgs doublet and the singlet can be 
probed by the upcoming experiments at the LHC. 

The neutrino masses can be introduced by means of the 
following addition to the Standard Model lagrangian: 



C = C +N a (iP)N a -y aa HL a N a 
-^SN c a N a + V(H,S) + h.c, 



(1) 



where £0 includes the gauge and kinetic terms of the 
Standard Model, H is the Higgs doublet, S is the real 
singlet, L a (a — e,/i,r) are the lepton doublets, and N a 
(a = 1, ...,n) are the additional singlet neutrinos. Let us 



consider the following scalar potential: 

V(H,S) = m 2 |#| 2 +m 2 S 2 + A HS |iJ| 2 S 2 
+A H |ff| 4 . 



A C S 4 



(2) 



After the electroweak symmetry breaking, the Higgs 
doublet and singlet fields each develop a VEV, (H) = 
vq = 247 GeV, (S) = V\, and the singlet neutrinos ac- 
quire Majorana masses M a = h a v±. The masses of the 
Higgs doublet and singlet at zero temperature are 



X H Vo+X s vf± 



(3) 

Below the scale 
of this symmetry breaking, the low-energy effective la- 
grangian is 



where D — (A h «q — A s w 2 ) 2 + A 



£ = £ SM + N a {t0)N a - 



- y aa H L a N a - ^ N c a N a + h.c. 



(4) 

where £sm is the Standard Model lagrangian. This is the 
usual seesaw lagrangian The number n of singlet 

neutrinos is not limited by the anomaly constraints or 
any other theoretical considerations, and the experimen- 
tal limits exist only for larger mixing angles [141 ]. Super- 
nova 1987A provides a constraint [7J, |l5|, which depends 



on the mixing angle and the sterile neutrino mass 
To explain the neutrino masses inferred from the atmo- 
spheric and solar neutrino experiments, n — 2 singlets 
are sufficient (l6j |. but a greater number is called for if 
the lagrangian ((4]) is to explain the LSND [lTfl. the r- 
process nucleosynthesis [18|, the pulsar kicks [2j, 
and dark matter 0, 0, BE Gl • The scale of the right- 
handed Majorana masses M a is unknown; it can be much 
greater than the electroweak scale [l3| . or it may be as 



low as 1-10 eV [17|, [20(. The seesaw mechanism [13j can 
explain the smallness of the neutrino masses in the pres- 
ence of the Yukawa couplings of order one if the Majo- 
rana masses are much larger than the electroweak scale; 
then the light neutrino masses are suppressed by the ratio 
(H)/M. However, the origin of the Yukawa couplings re- 
mains unknown, and there is no evidence that these cou- 
plings must be of order 1. In fact, the Yukawa couplings 
of the charged leptons are much smaller than 1. Theo- 
retical naturalness arguments in favor of the low-energy 
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seesaw [l?} appear to be as compelling as those in favor of 
the high-scale seesaw [ljj. In both limits one can have a 
successful leptogenesis: in the case of the high-scale see- 
saw, the baryon asymmetry can be generated from the 
out-of-equilibrium decays of heavy neutrinos 21| , while 
in the case of the low-energy seesaw, the neutrino oscil- 
lations produce the asymmetry 22j. 

Let us first consider just one singlet neutrino, with 
Majorana mass of the order of several keV, in light of its 
possible role in explaining pulsar kicks and dark matter. 
The same neutrino can play an important role in star 
formation and reionization of the universe [B| , as well as 
other astrophysical phenomena [23| . 

The range of parameters consistent with the explana- 
tion of pulsar kicks is shown in Fig.[TJ It contains two re- 
gions corresponding to the resonance oscillations Q and 
off- resonance oscillations Q. The boundaries of these 
regions are defined by the requirements of (i) the super- 
nova energetics, including the SN1987A constraints, (ii) 
adiabaticity and weak damping for resonant oscillations, 
(iii) sufficient anisotropy to explain the observed pulsar 
speeds [ISl. 

If the coupling h in eq. (Q} is small enough, the sterile 
neutrinos are out of thermal equilibrium at any time after 
inflation. This is the case if the annihilations NN — > 
NN, NN — > scalars, etc. are not fast enough to keep the 
sterile neutrinos in equilibrium for temperatures T in the 
range rh s < T < T^eat- F° r rh ~ 1 TeV and the reheat 
temperature T re heat < 10 16 GeV, the sterile neutrinos are 
out of equilibrium for h <C 10 -6 [121 ]. 

The sterile neutrinos can be produced out of equilib- 
rium in two different processes. First, since S is in ther- 
mal equilibrium at high temperature, sterile neutrinos 
are produced through decays S — ► NN ll|. Most of 



these neutrinos are produced at temperatures of the or- 
der of rh s ~ (0.1 — 1) TeV. Second, at much lower temper- 
atures, T ~ 0.1 GeV, the sterile neutrinos are produced 
from oscillations of active neutrinos, as in the Dodelson- 
Widrow (DW) scenario If the lepton asymmetry is 
relatively large, the sterile neutrinos are produced much 
more efficiently |9j. The lepton asymmetry of the uni- 
verse is not know, but strong upper bounds do exist 24 1 . 

These two production mechanisms operate sequen- 
tially. S decays are governed by the coupling h, while 
the production via DW mechanism depends on mixing 
angles. In the limit y a iVo <C min^ (Mi), these mixing 
angles are given by the usual seesaw relations: 9 a i = 
Uai(H) /Mi. For simplicity we will assume that only one 
sterile neutrino has mass of several keV and that only 
one of the mixing angles is non-zero. This mixing an- 
gle is a function of both couplings, y and h, as well as 
tan/3 = wq/ u i : 




(5) 



There are two limits which, in combination, appear to 



FIG. 1: The x-ray limits reported in Ref.[13] (dashed line) 
apply if the sterile neutrinos account for all the dark mat- 
ter (yi s = 0.2). The value of Sl 3 depends on the production 
mechanisms, but it cannot be lower than the amount pro- 
duced via DW mechanism 0) (except for the low-reheat sce- 
narios |2p|. [29l] ) . The model-independent exclusion plot (pur- 
ple region) is obtained by assuming this minimal value. A 
sterile neutrino with mass 3 keV and sin 2 6 ~ 3 x 1CP 9 , pro- 
duced at some temperature above 100 GeV, can explain both 
pulsar kicks and dark matter. 



stymie the DW proposal 0] for sterile dark matter. We 
will see that the models with singlet Higgs bosons offer a 
way to reconcile these seemingly contradictory bounds. 

One important limit on the abundance of the relic ster- 
ile neutrinos comes from the x-ray observations. The relic 
neutrinos can decay into one of the active neutrinos and 
a photon. Since this is a two-body decay, the photon en- 
ergy is equal one-half of the particle mass. The flux of 
expected x-ray photons depends on the decay rate, which 
is proportional to (sin 2 6) [25j . This flux also depends on 
the relic abundance of sterile neutrinos, which, in turn, 
depends on the value of h (for the S — ► NN production 
at temperature T ~ m s ) and 6 (for the production via 
neutrino mixing at temperature T ~ 0.1 GeV). If the 
latter is the dominant production mechanism, then the 
sterile neutrinos can be the dark matter in the range of 
masses [26] below 3.5 keV, as shown in Fig. [T] For larger 
masses, the x-ray limits on the decay photons 27J disal- 
low the relic sterile neutrinos, unless their abundance is 
smaller than what is required for dark matter. However, 
if S decays at T~ (0.1 — 1) TeV are the dominant source 
of sterile neutrinos, they can still make up the entire dark 
matter, even if 9 is arbitrarily small. 

Another limit on dark matter comes from the Lyman- 
a observations (28j and is based on the requirement that 
the dark matter be cold enough to generate the smallest 



3 



structures observed in the absorption spectra of distant 
quasars. This limits the free-streaming length of dark 
matter, whose relation with the particle mass depends 
on the production mechanism. If the sterile neutrinos 
are produced via mixing, then the Lyman-a bound is 
m > 10 keV (28[. Obviously, this limit is in conflict with 
the x-ray limit mentioned above, namely m < 3.5 keV. 

If the sterile neutrinos account for only a part of dark 
matter, then the Lyman-a bounds do not apply, and the 
x-ray bounds are weaker. The x-ray limits 27J, shown 
as a dashed line in Fig. [1] are based on the assumption 
that the sterile neutrinos constitute the entire dark mat- 
ter, i.e., £l s = 0.2. However, if h — ► 0, this condition 
is not satisfied for points along the dashed line; in fact, 
il s -C 0.2 for most of these points. This is because, in 
the absence of S decays, the density of sterile neutrinos 
in the universe may be too small to explain dark mat- 
ter. The x-ray signal is correspondingly smaller in this 
case, and the x-ray bounds are weaker. However, there 
still exists an x-ray bound on a sterile neutrino with a 
given mass and mixing angle, although this particle may 
be only a small part of dark matter. 

The production via neutrino mixing [|| provides the 
lowest possible abundance of sterile neutrinos, as long as 
the universe was reheated to a GeV or higher tempera- 
ture. In low-reheat cosmological scenarios the bounds are 
relaxed considerably 2(1 29|. Here we will assume that 
the universe has reached temperatures well above TeV 
after inflation. In Fig. Q]we show both the bounds based 
on the assumption fl s = 0.2 and the model-independent 
exclusion region (solid purple) based on the production 
only via the DW mechanism [6j, in which case fl s can be 
smaller than 0.2. To calculate this production, we used 
an analytical fit to the numerical calculation of Abaza- 
jian [261 ] . There may be some hadronic uncertainties in 
this calculation [3Cj |. 

There is a range of parameters for which the sterile 
neutrinos can explain the pulsar velocities and can af- 
fect the star formation, although they may not be the 
dominant component of the dark matter (see Fig. [IJ. 
However, it is also possible to explain pulsar kicks and 
the dark matter simultaneously. Indeed, if the S'-decay 
mechanism dominates the production of relic neutrinos, 
they may be redshifted sufficiently as the universe cools 
from T s - (0.1 - 1) TeV to T < MeV. 

If all the dark matter is made up of sterile neutrinos, 
the bound m s > 10 keV 28] applies to the DW scenario, 



in which the average momentum of a keV neutrino at low 
temperature T is (p s )i)W ~ 2.8 T [2f|. Any additional 
redshifting of sterile neutrinos (e.g., due to some entropy 
production [^H) relaxes the 10 keV limit to a lower value. 

Production of sterile neutrinos via decay S — ► NN 
occurs mainly at temperature T s w fn s [ill ]. Our model 
has two mass eigenstates given by eq. ((3} , subject to ther- 
mal corrections. In the case of different masses and non- 
negligible mixing, one of the mass eigenstates decouples 



before the other. Each mass eigenstate has a non-zero S 
component and, therefore, contributes to the production 
of the relic population of sterile neutrinos. If we make a 
simplifying assumption of equal masses (achieved when 
X U (H) 2 i=s X S {S} 2 and when A HS is small), the results of 



Ref. [ll| are directly applicable to our case, and one can 
write that the amount of sterile dark matter as 



fi., = 0.2 



•<s> 



33 



£ / \1.4 x 10~° / \ ,// . 



, . 33 
0.2 — 



£ J \lAx 10- 



(6) 
(7) 



where £ is the change in the number density of sterile 
neutrinos relative to T 3 due to the dilution taking place 
as the universe cools down from T s ~ (0.1 — 1) TeV 
to a temperature well below MeV. The sterile neutrinos 
produced in these decays have an almost thermal spec- 
trum at the time of production. More precisely, their 
average momentum (p s ) = tt 6 /(378C(5))T « 2.45 T [n|], 
while the average momentum of the relativistic fermions 
in equilibrium is p T ~ 3.15 T. As the universe cools 
down, the number of effective degrees of freedom de- 
creases from g*(T s ) = 110.5 to g*(0.1MeV) = 3.36. 
Then £ = g*(T g )/0*(O.lMeV) « 33. This causes the 
redshifting of (p s ) by the factor ^Z 3 : 



(P 



si(T«lMcV) 



0.76 T 



110.5 



g*(m s ) 



1/3 



(8) 



Comparing eq. {8]) with the DW case, one concludes that, 
as long as the population of sterile neutrinos is dominated 
by those produced at high temperature (large enough h, 
small 8), the Lyman-a limit changes from 10 keV to 



> 2.7 keV 



(9) 



This lower bound is shown in Fig. [TJ 

In addition, the electroweak phase transition in this 
model can be first-order, and the entropy production 
can further redshift the momentum distribution of sterile 
neutrinos. Enqvist et al. [12] have found that the energy 
density increase due to the phase transition can be at 
most 10 Tq, where To is the transition temperature. This 
changes £ by at most a factor 1.3. 

The presence of the singlet Higgs boson has important 
implications for the Higgs searches. Although we have fo- 
cused so far on the keV sterile neutrino, the present data 
requires n > 3 sterile neutrinos in eq. dU) [13) El • This is 
not in conflict with big-bang nucleosynthesis (BBN) : one 
or two additional thermalized neutrinos are consistent 
with the BBN constraints at one (two) sigma level [32j |: 
besides some of these sterile neutrinos can be out of equi- 
librium. The couplings h a of the additional sterile neu- 
trinos need not be small, and the mixing between H and 
S can also be large. Hence, the Higgs boson can decay 
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invisibly. This possibility has the effect of weakening the 
LEP bound on the mass of the lightest Higgs [H, [3^ . 
However, one can discover the invisible Higgs /i; nv at the 
LHC in the Z + h- lm channel 34j. For Higgs mass of 
120 GeV, the discovery is possible at the LHC already 
with 10 ftr 1 in the Z + h im channel, while 100 fb" 1 of 
data can afford the discovery in the weak boson fusion 
channel [34|. For some range of couplings, the singlet S 
can decay into the visible channels and can be discovered 
via displaced vertices [HI]. 

To summarize, the relic sterile neutrinos with mass 
of a few keV can simultaneously explain the pulsar ve- 
locities and dark matter if their production in the early 
universe occured above the electroweak scale. Over a 
broader range of parameters (Fig. [1]), the sterile neutri- 
nos can explain pulsar kicks and can play a role in the 
star formation, while they may not be the dominant com- 
ponent of dark matter. 
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